Bosonization for Wigner-Jordan-like Transformation : 
Backscattering and Umklapp-processes on Fictitious Lattice. 
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We analyze the asymptotic behavior of the exponential form in the fermionic density operators 
as the function of ruling parameter Q. In the particular case Q = n this exponential associates 
with the Wigner- Jordan transformation for XY spin chain model. We compare the bosonization 
approach and the evaluation via the Toeplitz determinant. The use of Szego-Kac theorem suggests 
that at Q > ir/3 the divergent series for intrinsic logarithm provides a bosonized solution and faster 
decaying one, found as the logarithm's value on another sheet of the complex plane. The second 
solution is revealed as umklapp-process on the fictitious lattice while originates from backscattering 
terms in bosonized density. Our finding preserves in a wide range of fermion filling ratios. 
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The method of bosonization proved to be an effective 
tool for studying the properties of one-dimensional (ID) 
systems of interacting fermions. Describing the fermion 
system through the bosonic collective modes, 1 this tech- 
nique was initially used for providing the exact solu- 
tion of Luttinger model. 2 Recently the bosonization at- 
tracted a renewed interest due to investigation of Quan- 
tum Hall Effect 3 , two channel Kondo model 4 and rIigh-T c 
superconductivity. 5 Essential efforts are made to gener- 
alize the method for higher dimensions. 6 

In a recent paper 7 the bosonization technique was 
successfully employed to study the magnetic properties 
of the high-T c -related system. The ID model with a 
spin-charge separation was solved and important physi- 
cal consequences were obtained. The bosonization was 
used to estimate the (time-dependent) Wigner-Jordan- 
like exponential in fermionic on-site density operators 
CjCj. Namely the quantity of the interest was the av- 
erage ^exp iQY^j=i c ] c j) a ^ arrj itrary value of Q. This 
average should be 2ir— periodic in Q and attain the real 
value at Q — ir (the case for Wigner-Jordan transfor- 
mation). These properties are lost after the simple use 
of the bosonization technique. Although various ad hoc 
tricks were implemented to cure the problem at Q = n 
(see, e.g. 8 and references therein), no systematic bosonic 
representation existed. 

In an attempt to resolve this puzzle we addressed our- 
self to the vast available literature on the subject. The 
results of this survey could be outlined in the following 
way. The straightforward application of bosonization, 
though formally valid, is seemingly incapable to provide 
a closed form of the final answer due to the complexity of 
intermediate expressions. The key could be found in com- 
parison of two counterparts for the same problem : the 
bosonization by Mattis and Lieb 2 and Luttinger's origi- 
nal solution 9 via the Toeplitz determinants. The Szego- 
Kac theorem, used in the approach by Luttinger, dates 
back 10 to Onsager famous solution of 2D Ising model. 
The exhaustive literature on the latter subject suggests 



some useful hints about the discussed issue. 

Basing on the alternative representation of the consid- 
ered average, we make a following major statement. The 
above exponential form in the fermionic density opera- 
tor reveals a singularity at some finite value of Q. This is 
verified by an application of Szego-Kac theorem and man- 
ifests itself as the appearance of umklapp-processes con- 
tribution on the fictitious lattice (see below). We provide 
an explicit expression obeying the above 27r-periodicity 
and real valued at Q — w. 

The rest of the paper is organized as follows. We for- 
mulate the simplified version of our problem and solve 
it by bosonization method in Section I. We rewrite our 
quantity of interest as the Toeplitz determinant and ap- 
ply the Szego-Kac theorem in Section II ; the appearing 
singularity is discussed here. In Section III the structure 
of the desired solution is suggested and justified. This 
solution and its physical consequences are discussed in 
the Section IV of the paper. 



I. THE MODEL AND APPROACH BY 
BOSONIZATION 

We consider a gas of free spinless fermions on a lattice 
with Hamiltonian 



Hi = ^^(cosfci? — cos p)c|c 



(1) 



where Fermi wave- vector kp = en. Our aim is to evaluate 
a following quantity : 



C(Q,l) = (exp iQJ2 c ] l 



(2) 



Hi 



where CjCj corresponds to the on-site density operator. 



The expression of the form (2) appeared previously in 
various studies concerning the physics in one dimension. 
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For the particular case of Q = tt it is found in the Wigncr- 
Jordan transformation for the spin operators in ID XY 
model. It can be also associated with the incorpora- 
tion of interaction effects in the bosonization approach 
to Luttinger model. In our problem (see Appendix A for 
details) the quantity (2) connected the observable spin 
susceptibility and easier tractable one on the fictitious 
lattice. The parameter Q appeared as a wave-vector be- 
longing to the fictitious lattice and varied over the whole 
Brillouin zone, Q € [—tt,tt). 

The average (2) is easily evaluated by introducing the 
boson representation for the fermionic density. The main 
steps of this derivation are outlined below. First we in- 
troduce the density operator : 



c l> c p+q> 



where N is the total number of lattice sites. 



JQcl 



exp 



(3) 



(4) 



and we have defined 4>(l,q) = (1 — e tql )/(l — e lq ) and 
<j>(l,0) = I. The non- fluctuating component of the density 

(cjCj) —cis treated separately and leads to the prefactor 
e iQd m g q 4 

Now we may use a well known procedure 11 in which 
the fermion density is split into contributions from right 
and left moving fermions. We write 

p{q) = p+(q) + p-{i) 



p>0 



p(q) = N- 1 (N_+N + ), 



p<0 
q = 



(6) 



where N + (A_) are charge operators for right (left) 
movers, describing the fluctuations of the Fermi surface 12 
The long wave-length density fluctuations of the I — 
d electron gas may be represented by boson operators 
defined by the relations : 



p+(q) = v(q) {6(q)b\ + 6{-q)b ) 
p-(q) = v(q) (6( q )b_ q + 9(- q )bl), 



(7) 



where v(q) = y/\q\/(2Nw). 2 In terms of the above oper- 
ators, the Hamiltonian reads : 



H = v F ^\q\blb q + v F ^(Nl+Nl) 



(8) 



where Fermi velocity v F = s'mk F . 

With these definitions we may rewrite Eq. 4 in the form 



C{Q,l)=e i ^ l (e^i^{N + + N_)^ 

x(e X piQj2Hq)^q,l)[bl + b_ q ]} (9) 



Let us focus on the case of zero temperature; the finite T 
are considered in the Appendix B. The fermionic vacuum 
corresponds then to the absence of bosons in the ground 
state, in addition N± = 0. The evaluation of the above 
expression yields immediately : 



C(Q,l) =exp 



iQcl 



91 

2tt 2 



Hi) 



where 



a 



1 — cos q(l — I') 
1 — cos q 



(10) 



(11) 



The integration in Eq. 11 is cutoff at a wavevector of 
order k F since the density fluctuations of the electron gas 
are well defined excitations only in the region \q\ < k F . 13 
Evaluating the above integral in the logarithmic approx- 
imation we find : 



C(Q, I) OC CbosoniQ, I) = CXp 



Q 2 

iQcl 7T In I 

2tt a 



(12) 



This expression indicates the power law decay of the cor- 
relations with I, which is well established fact for the ID 
systems. 14 

The result (12) seems however unsatisfactory on two 
following reasons. First, the initial expression (2) was 
invariant under the shift Q — > Q + 2n. This property is 
obviously absent in (12). Second, e 47rc 3 Cj = ±1 for each 
j, and one should have the real value of (2) at Q = n, 
meanwhile this feature is also not found in (12). 

The inadequacy of the expression (12) for the partic- 
ular case Q — tt was noticed by Luther and Peschel 11 . 
They argued that the term oc e ^^ cl should be added 
to (12) to get the correct result. The necessity of inclu- 
sion of such term was ascribed to the 2k F ("backscat- 
tering") processes ( we remind that 2k F — 2nc and 
Qc — 2k F = —Qc). Really, one might argue that the 
representation (5) of the electronic density is incomplete 
and the low-energy fluctuations with q ~ 2k F should 
be explicitly included into the notation. 12 These terms 
are described in the Appendix C where we show that the 
"2k F " terms produce very complicated and seemingly in- 
tractable expressions. 

A general question arises here. If 2k F processes pro- 
vide the real value of C(Q, I) at Q — tt, then why these 
2k F terms are less important at other values of indepen- 
dent parameter Q where C(Q,l) is not real ? Below we 
argue that at Q > tt/3 the "marginal" contribution of 
backscattering breaks the analiticity of C(Q, I). 

In our problem the parameter Q varies from to tt. 
Therefore the ignorance of the above discrepancy might 
produce non-trivial consequences. To resolve this ques- 
tion we formulate below the alternative way for the eval- 
uation of C(Q,l) (2) and provide the answer which ex- 
plicitly obeys two above conditions, the 27r-periodicity in 
Q and real valued quantity C(Q = tt, I). 
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II. SZEGO-KAC THEOREM 

We use a modification of an approach originally due to 
Lieb, Schultz and Mattis 15 . It is based on a representa- 
tion of the correlation function C(Q, I) of Eq. (4) as the 
Tocplitz determinant 



C(Q,l)=det(M^) 



(13) 



5 km - (1 - e lQ )g, 



km • 



The I x I matrix is constructed out the fermionic 
correlation functions 



9kn 



( c l c J - 2n 



sm(nc(k — m)) 
Tr(k — m) 



dqn q e lq{k - m) 



(14) 



where the Fermi function n q = 1 at \q\ < fop and n q = 
otherwise. 

For completeness we sketch below the derivation of 
(13). We observe that 

e *0^ = ( c t+ c .)( c t +e «J c .) (15) 

which is readily verified in the representation diagonaliz- 
ing c\ci 



• ■ Defining 



A l = c\+c i , Bi = c} + e 



■t 



we have 



C(Q,l) = (A 1 B 1 A 2 ...B l _ 1 A l B l ) 



(16) 



(17) 



Zuber and Itzykson evaluated the square of the value (2) 
in a similar problem. 16 They noticed at Q — ir one had 
the possibility to identify the operators Ai, Bi with a new 
fermionic (Majorana) field. This is not so for other values 



of Q since AiBi + BiAi = 1 



JQ 



^ 0, therefore the trick 



used in 16 with introducing the second Majorana field is 
hardly applicable to our case. 

However one can apply the Wick theorem which per- 
mits one to express the vacuum expectation value (17) 
in terms of expectation values of products of just two 
operators. 17 This procedure is facilitated if one notices 
that 



(AiAj) = (BiBj) = 0, i ? j 



and 



(1 



e iQ )g l3 



(18) 



(19) 



The most straightforward contribution to (17) is 

(A 1 B 1 )(A 2 B 2 )...(A l B l ) 

All other pairings are obtained by permuting the B's 
among themselves with A's fixed. 18 The number of cross- 
ings of ^4's by B's is always even, hence the sign associ- 
ated with a given permutation is (— l) p where p is the 
signature of the permutation P of the B's. As a result 



C(Q,l) = ^(-l) p ' ' (A 1 B P[1) )(A 2 B P(2) ) . . . (A t B P{l) ) (20) 
p 

and we arrive to the above Eq. (13). 

Note that the representation (13) docs not appeal to 
any particular law of fermionic dispersion, and the notion 
of the Fermi wave- vector only appears in g^ m (at T = 0) . 

Now we use the asymptotic behavior of the above de- 
terminant at large I known as Szego-Kac theorem 19 ' 10 , 



lim 



where we follow the notation by Luttinger 9 
D = cxp / 2^1n/„ 



Ki 



r *± 



'In /a, 



Where 



f q = J2 J q(m - k) Mi2\ 



(21) 

(22) 
(23) 

(24) 



m— — oo 



in our case f q = 1 — (1 — e 



Therefore at the first 



glance one has ln/ 9 = iQn q . Next, D = iQc, Ki = iQgio 
and simple calculation returns us back to Eq. (12). It was 
noticed in 9 however that the proof of Eq.(21) required the 
convergence of the series for logarithm 



w 9 = -E 



(!-/«)" 



n=l 



n 



(25) 



In his discussion Luttinger explicitly assumed the conver- 
gence of the above series, associating it with the smallness 
of the interaction between electrons 9 (see also 2 ). This ar- 
gument is not applicable to our case since 



ln/g = -n q J2 



71=1 



(-2isin(Q/2))" c<Qn/2 
n 



(26) 



Although one might say that this series converges to the 
value iQn q , it is divergent in absolute values at \Q\ > 
it/3. Thus formally it can produce an arbitrary value as 
its sum upon the reordering of terms. The consequences 
of this observation are elucidated below. 

Concluding this section we wish to make some remarks 
concerning Eq.(21). 

The attentive reader may notice, that the proof of 
Szego-Kac theorem 19,20 usually demands that the func- 
tion f q is continuous. In the absence of interaction this 
is achieved at finite temperatures, when the Fermi distri- 
bution function is smeared at the scale q~ kp ~ T/vp. 
In this latter case however one expects the exponential 
tail 7 as a far asymptote of Cq(1) and power-law decay 
as intermediate regime. The problem is hence somewhat 
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similar to the problem of correlations in the scaling limit 
(uniform asymptote) in the 2D Ising model near T c . 21 On 
the other hand, the sufficient conditions for the validity 
of Szego-Kac theorem are apparently not known. 20 The 
scrupulous investigation of the arising difficulties would 
bring us far from our initial task. That is why we tried to 
get an insight of the situation by the procedure described 
in the next section. 



III. THE CONJECTURED SOLUTION 

We numerically calculated the determinants of I x I 
matrices of the form (13) for different values of c and Q 
and for 1 < I < L = 64. We found that the asymptotic 
behavior det(M^) oc l~® l 2lx actually begins with I of 
order of unity. We illustrate this fact on the Fig. 1 where 
we plot the auxiliary quantity 



A(Q,l) =det{M^)l Q2 ^ 2 



(27) 



for some particular values of Q and c. If Eq. (12) would 
hold, | A(Q, I) | would be constant, which is indeed realized 
at largest I. An essential additional feature could be also 
noted on this figure. Namely the extra oscillating term 
exists, which decays faster than the bosonized solution 
(12) at small Q. However with the increase of Q, the 
decay of this extra term becomes comparable to one of 
(12) and both terms are of the same significance at Q = 
it. This qualitative picture is reproduced at the other 
values of c as well. 

It is instructive to analyze at this step the Fourier spec- 
trum of the above amplitude : 



A(Q,k) = J2e~ iklA (Q,l) 



(28) 



Again if our initial guess about the dependence (12) in 
the whole range of Q were true, we would observe the 
following simple behavior 



A^(Q,k)oc5(k-cQ) 



(29) 



One can see on the contour map of A(Q, k) (Fig. 2) that 
at small Q our initial expectations are satisfied. At the 
same time at larger Q the second solution comes into play 



A( 2 \Q,k)ot5(k-c(Q-2ir)) 



(30) 



This finding shed some light on the nature of the dis- 
cussed discrepancy. Really the "umklapp" term A^ may 
take its origin in the multi-sheet structure of the above 
logarithm (26) on the complex plane. 

At Q > 1 when the series (26) diverges, the value of 
ln/ g = n 9 lne 4< 2 could be both iQn q and i(Q — 2Tr)n q , 
i.e. the value of logarithm on the other sheet of complex 
plane. 

We conjectured that the correct answer may be a lin- 
ear combination of the bosonized solutions Cb OS on(Q,l) 



(12) and Cb OS on(Q — 2tt, I). We tried to fit therefore the 
I— dependence of C(Q, I) at fixed c and Q : 

C{Q, I) fit = ai(Q)C boson (Q , I) + a 2 (Q)C boson (Q - 2tt, I) 

(31) 

with the real functions a.\ 2 {Q) as the parameters of the 
fit. 

The results of this fit are encouraging. The dependence 
(31) and the calculated values (13) are indistinguishable 
at I > 1, < Q < 7T and at 0.2 < c < 0.8. 

The above functions a\ y2 {Q) slightly vary with c in the 
above interval, and are symmetrical with respect to the 
point c = 0.5 as one could expect in view of the particle- 
hole symmetry of the initial Hamiltonian. The typical 
variance in Q is shown on the Fig. 3a for c = 0.5. 

One sees that ot\{Q = 0) = 1 and a 2 (Q = 0) = 0, while 
the ending values cci(7r) and 02 M coincide so as to pro- 
duce the real values of C(Q = ir,l). We noticed next 
that the slopes of ai^(Q) at Q = it coincide differing in 
sign. We redraw the obtained picture by shifting it by 
±2*7r as shown on Fig. 3b. Recalling the 27r-periodicity 
of the initial expression (2) it becomes evident that actu- 
ally a±(Q) and a 2 (Q) are the same even function a(Q), 
obeying the relation : 



a(Q) = ai(Q), Q < ir 

= a 2 (2ir-Q), Q > tt 



(32) 



We now come to our main hypothesis. In view of the 
multisheet structure of the logarithm found in the expres- 
sion for the asymptotic behavior of Toeplitz determinant, 
we conjecture that the asymptotic behavior of the aver- 
age (2) can be represented in the following form 

00 

C(Q,l)= C boson (Q-27rm,l)a{Q-2irm) (33) 

m— — oc 

We notice that, numerically, the function a(Q) is close 
to Gaussian, 

a{Q) ~ cxp ^2 ln/ , 

with Iq < 1. It allows one to further compact the expres- 
sion (33). Using the definition of elliptic theta- function 
$3(2:, q) 22 , we write an approximate equality (cf. Eq.(12)) 



C(Q,Z)~exp 



Q 2 I 
iQcl - — 2 In — 
2-k 1 lo 



iQ I I 2 
ircl H In-, -r 

TT Iq l Z 



(34) 



which holds for the arbitrary Q. 

Let us discuss here the particular consequences of 
Eq.(33) in view of our previous findings. 7 . Using the 
bosonization approach we determined the dynamics of 
the fermionic correlations on the partially filled chain; 
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this result could hardly be obtained by another tech- 
nique. Although we were provided only by the bosonized 
solution of the type (12), our main conclusions in 7 pre- 
serve as explained below. 

Let us write the bosonized solution in the form 
Cboson{Q,l) — e z ® cl F(Q,l), the Fourier transform 
F(Q,k) ~ \k\ Q2 / 2v2 - 1 . Then the quantity entering the 
formula (A4) takes the form 

C(Q, Q-k) = J2 [ ^ F( - Q ~ 2lTm > ~ 2rrm ) 

m,n 

x6[(Q-2nm)(l + c)-k-2im + q]. (35) 

The terms in (33) with m ^ stemming presumably from 
backscattering in real space are manifested in (35) as the 
umklapp processes on the fictitious lattice. They restore 
the invariance of the observable quantity x(fc) (A4) upon 
the shift Q — > Q + 2tt, i.e. upon the choice of the Brillouin 
zone on the fictitious lattice. This corresponds to the 
assumed cyclic boundary conditions. 

It was shown in 7 that x(&) had a principal contribution 
from Q ~ ±7r, q ~ 0. From the form of a(Q) and the 
restriction \Q\ < it it is clear that only the terms m = 
0,±1 are present in (35). A simple analysis shows then 
that the only modification of our previous results would 
be a symmetrization of the left and right shoulders of the 
function \{k) (A4) around its peak values, (see Fig. 3 in 
Ref. 7 ) 

IV. DISCUSSION 

Firstly one notices that (33) is explicitly 2ir periodic 
in Q. In view of evenness of a(Q) the Eq. (33) is real at 
Q = 7r as it must. 

Secondly, the leading asymptotic behavior at \Q\ < ir 
is obviously delivered by m — term in (33), which is 
the term produced by the bosonization procedure. At the 
same time, at \Q\ > 1 the divergence of the underlying 
logarithm series (26) manifests itself in the appearance 
of the "neighboring" in Q solution. This second solution 
is relatively unimportant at 

cxp(-^ — InZj <exp(-l) 

therefore one can ignore it only for exponentially large 
distances 

Z>exp[7r/2(7r-Q)]. (36) 

This argument can be also applied to the Fig. 2 where 
the second component is present mostly at tt — Q < 
7r/(2 In I,) ~ 0.4 

Thirdly, the bosonization is usually expected to work 
well, if the filling ratio c is close but not equal to 0.5. 
This is to preserve the almost linear character of disper- 
sion near the Fermi level and to avoid the influence of 



2kp processes. We saw that the Eqs. (33), (34) held 
for a rather wide range of c, hence 2k f processes were 
important for all c, while the bosonized solution was a 
robust core in the final expressions. 

Fourth, the algebra of the operators b q , 6j, N±(U±) is 
complete 12 and a proper treatment should produce con- 
sistent expressions at any Q. This could be an unfeasible 
task in view of the complexity of expressions sketched in 
Appendix C. Here we note that since (33) represents the 
asymptotic behavior, a(Q) could be non-analytic func- 
tion. This point and the very structure of Eq. (33) im- 
plies that the thorough consideration of backscattering 
terms in the fermionic density operator can lead to (33) 
only by careful resumming of the intrinsic divergences in 
intermediate expressions. 

Thus we observe that Eq. (33) does not contradict to 
the previously known results while its detailed analytical 
verification may demand enormous efforts. 

At the moment it is unclear whether our result should 
be regarded as a very particular one or it could be found 
in other models, too. The power-law decay of correla- 
tions is usually analyzed in terms of leading asymptote, 
the existence of faster decaying one, as it is in our case, 
may be missed. This should be compared to the results 
of Ref. 23 , where the faster decaying asymptotes of corre- 
lations were analyzed. It was shown, in particular, that 
these next-to-leading terms may dominate for the models 
with finite-range fermionic interactions. 

In conclusion we show that the Wigncr- Jordan-like ex- 
ponential form in fermionic density operators reveals sin- 
gularity at some finite value of ruling parameter. The 
appearing additional terms, apparently stemming from 
the 2/cf processes on real lattice, correspond to umklapp 
processes on the fictitious lattice. These terms restore a 
certain symmetry of the problem, lost after the simple 
use of the bosonization technique. 
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APPENDIX A: THE CONNECTION BETWEEN 
FICTITIOUS AND REAL LATTICES 

Let us consider a chain with 2N sites and N(l + c) 
electrons on it, we assume N ^> 1, < c < 1. The sites 
labeled by integer index are always occupied by electrons. 
The remaining cN electrons reside on N sites with half- 
integer indices. No double occupancy is allowed. The 
elementary kinetic process (disregarding the spins) is the 
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movement of electrons on the half-integer sites. This pro- 
cess can be thought as combined of two steps, first one is 
to allow the electron from the site i to hop onto the neigh- 
boring (empty) half-integer site i + 1/2. At the second 
step an electron moves from i — 1/2 to i. One sees that 
upon this processs the spin sequence of the whole chain 
preserves. One can also adopt that the magnetic ex- 
change between all neighboring spins has the same value. 

For such a system, it is convenient to introduce the 
notion of fictitious lattice of spins with (l+c)N sites. The 
spin dynamics on the fictitious lattice is then described 
by the Heisenberg Hamiltonian 



N+cN-1 

H 2 = J ^2 a-i.er, 

i=l 



(Al) 



while the kinetics of the system is that of the fermionic 
gas (1). 

Despite of the simple form of H = Hi + Hi the ob- 
servable spin-spin correlation function, 



X (k) = (S fc .S_ fe ) 



(A2) 



is non-trivial because the relationship between Si and er^ 
is not simple 24 . The difference between the coordinates 
of a integer-labeled spin measured in the real and in the 
fictitious lattices is given by the number of fermions lo- 
cated to its left. In terms of er fc , the Fourier transform of 
the spins on the fictitious lattice, we have: 



1 



VN + cN 



^2(tq cxp iQ(l + ^2^ j+1/2 c j+1/2 ), (A3) 



Substituting Eq. A3 in Eq. A2 and using the fact that the 
ground-state wavefunction of the decoupled Hamiltonian 
factorizes, x(k) may be written as a convolution 25 : 



i: 



dQ 
2tt 



f(Q)C(Q,Q-k), 



(A4) 



The first factor in Eq. A4, f(Q), is the spin correlation 
function of the fictitious ID Heisenberg model, 26 



f(Q) = (<7Q-<t-q)h 2 



(A5) 



The second factor, C(Q, k), contains the kinematic effects 
and is given by the Fourier transform of (2), where one 
can drop the factor 1/2 labeling spinless fermion position. 
Thus we see that the observable spin correlations x(k) on 
the real lattice are connected to those on the fictitious lat- 
tice f(Q) via the integral transform (A4) with the kernel 
C(Q,k). The principal contribution to (A4) is delivered 
by Q ~ ±7r where f(Q) logarithmically diverges 7 . 



APPENDIX B: TEMPERATURE EFFECTS 

At finite temperature we calculate the average of the 



form (9) via the basic identity : 



(A) = Tr (e- fiH A) /Tr (e~^ H ) 



(Bl) 



The eigenstates of the bosonized Hamiltonian are given 
by 12 : 



\N + ,N-,{n g }) = (U + r+(U-) 



(n,!)V2 



|0> 



(B2) 



Here U+ (U-) is the ladder operator which raises the 
fermion charge N + (iV_) in unit steps and commutes 
with bosons b q . 

It is then straightforward to show that the 
temperature-dependent correlator is given by the expres- 
sion : 



C{Q,1) = C{Q,1)t= F{(3 = 1/T,1) 



(B3) 



with 

FX/3, Z)= cxp 



exp(/3\q\v F ) - 1 



En exp[-(/37Tf F n 2 - iQln)/N] 
J2 n exp[-(3nv F n 2 /N] 



(B4) 



It is convenient to introduce here the temperature cor- 
relation length, £ = /3v F . The first factor in (B4) stems 
from the boson contribution and takes the form 



Q 2 sinh7rJ/£ 
exp J" In —— 



exp - 



in accordance with the result in 7 found with the use of 
Szego-Kac theorem. 

The second factor in (B4) stems from the charge fluctu- 
ations and is relatively unimportant for all temperatures. 
If £ 3> N then this second factor is unity and the temper- 
ature does not break the coherence within the length of 
a chain (or, more precisely, on a ring in view of implicitly 
assumed periodic boundary conditions). In the opposite 
case £ <C N we have the Gaussian decay of correlations 
: exp[— Q 2 ^ 2 /(47r^7V)] The appearing new length scale 
\/£,N 3> £ is however irrelevant, since the correlations 
beyond £ are already damped due to the boson contribu- 
tion. 



APPENDIX C: "BACKSCATTERING" TERMS 

The consideration of the low-energy fluctuations of the 
electronic density p(q) with q ~ 2k f produces the extra 
term to the expression (5) 

c\c l = Y,e iql [ P+ {q)+P-{<l)\+r{l) (CI) 
t(1) - (jijWu+U^j*'® + /i.e.) (C2) 
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with the "current" variable 12 

*j(l) = k F l + jj(N + + N_) (C3) 

It is seen that C(Q,l) contains now the exponential 
from the linear combination and the exponential in bose 
operators. Presently we did not succeed in evaluating of 
this type of expression. 
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FIG. 1. The quantity \A(Q,l)\ defined in Eq.(27) shown 
as a function of the chain fragment length I. The plotted 
values are for c = 0.3 and for Q/ir = 0.45, 0.6, 0.75 from top 
to bottom, respectively; the lines are connecting the points. 
Note the appearance of faster decaying asymptote at larger 
Q. 
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FIG. 2. The absolute value of A(Q,k) denned in (28) 
plotted for the filling ratios c = 0.3,0.5 and 0.7, from top to 
bottom, respectively. The contours are drawn at 0.03, 0.05, 
0.1, 0.3 and 0.5. The appearance of umklapp solution (30) is 
seen at large Q. 
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FIG. 3. a) The amplitudes of the solutions ai,2(Q) as 
determined by fitting the calculated quantity A(Q, I) with the 
Eq. (31) at c = 0.5. b) The same functions shown in the 
extended range of Q. 
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